Viscous shear-banding in foam 
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Shear banding is an important feature of flow in complex fluids. Essentially, shear bands refer to 
the coexistence of flowing and non-flowing regions in driven material. Understanding the possible 
sources of shear banding has important implications for a wide range of flow applications. In this 
regard, quasi-two dimensional flow offers a unique opportunity to study competing factors that result 
in shear bands. One proposal for interpretation and analysis is the competition between intrinsic 
dissipation and an external source of dissipation. In this paper, we report on the experimental 
observation of the transition between different classes of shear-bands that have been predicted to 
exist in cylindrical geometry as the result of this competition [R. J. Clancy, E. Janiaud, D. Weaire, 
and S. Hutzlet, Eur. J. Phys. E, 21, 123 (2006)]. 

PACS numbers: 83.80.Iz,83.50.Ax,83.60.Wc 



It has long been assumed that flowing and non-flowing 
regions can exist in driven complex fluids, especially in 
situations in which the stress is inhomogeneous. This 
behavior is often discussed in the context of yield stress 
materials [H, 0|. These are materials for which a crit- 
ical stress is required for flow to be generated. Nat- 
urally, in a situation in which the stress varies from 
values below to above the critical, or yield, stress, one 
will have coexistence of flow and no flow. More re- 
cently, experimental studies have revealed a rich va- 
riety of shear-banding situations that do not conform 
to the "traditional" yield stress picture of shear bands 
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One of the challenges to an increased understanding 
of shear banding in complex fluids is the large number 
of parameters and observed behavior. There is clearly 
no single cause of shear banding, and one must sort out 
the mechanisms for a range of different qualitative be- 
haviors. Among the interesting range of behavior and 
models, we are going to highlight a few cases because of 
their relevance to the work discussed here and because 
they highlight the importance of distinguishing between 
qualitatively different types of shear banding. It is im- 
portant to realize that these cases are not necessarily 
mutually exclusive. We will first consider issues related 
to the particulate nature of foams. 

When considering shear banding in many complex flu- 
ids, one has to determine if a continuum model will be 
relevant or if one has to consider in detail the fact that 
the fluid consists of macroscopic "particles" . For ex- 
ample, foam is gas bubbles in the micron to millimeter 
range surrounded by liquid walls. Granular matter is 
similar size particles with gas. Experimentally, a tran- 
sition to a "discrete" flow regime has been observed for 
extremel y na rrow shear bands on the order of 10 - 20 par- 
ticles fA Il7|. In this case, an interesting feature of the 
shear band is that the width becomes independent of the 
external driving velocity. This is fundamentally different 
behavior from continuum shear bands, where the width 



scales with the driving. 

Another qualitative signature of different types of 
shear localization is the continuity, or lack there of, of the 
rate of strain at the transition from flow to no-flow. Con- 
tinuous versus discontinuous transitions represent poten- 
tially different physics. This is supported by observations 
of a discontinuous transition in the rate of strain in a 
range of experimental systems [Til ] , including both three- 
dimensional systems and in quasi-two dimensional sys- 
tems (bubble rafts) @,E3]- For simple cases, the origin of 
the discontinuity may be understood in terms of a compe- 
tition between clustering of particles and the breakdown 
of cluster by shear [3] , but a complete theoretical under- 
standing of mechanism for selecting between continuous 
and discontinuous transitions is still lacking. 

Another qualitative feature of shear bands, especially 
in foams, is the fact that the velocity localization is of- 
ten correlated with the localization of other dynamical 
features in the system. For example, flow in foams is 
dominated by local rearrangements of the bubbles (Tl 
events). The degree of correlation between velocity and 
Tl localization, and the possible role of Tl events as a 
source of shear localization, are an important feature of 
shear bands. Models have suggested that shear bands can 
result when local stress inhomogeneities are long lived 
and produce spatio-temporal correlations in the occur- 
rence of plastic events in the flow (often associated with 
Tl events) 20]. This tends to result only in continuous 
shear bands, so it serves as an alternative model to ones 
that describe discontinuous transitions. More recently, 
a useful direction for exploring the role of Tl events in 
focusing the flow has been proposed using surface evolver 
[2ll ]. This work has shed interesting light on the role of 
polydispersity in controlling the degree of localization. 

As one starts to consider the various causes of shear 
localization, studies with quasi-two dimensional systems 
present an interesting option: role of external versus in- 
ternal dissipation. This case has been of interest for a 
class of systems involving a single layer of bubbles. For 
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these systems, there are three common configurations for 
the boundaries that are used to confine the bubbles to 
a plane. There is the classic bubble raft consisting of 
a single layer of bubbles at the air- water interface [22J. 
This has recently been labelled the liquid-air (LA) ge- 
ometry (23[ , which is a useful distinction when making 
comparison with other common quasi-two dimensional 
geometries for foams. One can also confine the bubbles 
between glass plates in a Hele-Shaw geometry, referred 
to as glass-glass (GG). Finally, one can study a mixed 
system using bubble on a water surface with a top glass 
plate, the liquid-glass geometry (LG). Detail studies of 
the dynamics of individual bubbles in these geometries 
have helped our understanding of the differences between 
the systems j23j. Experimental studies in all three sys- 
tems have established conclusively that dissipation be- 
tween the bubbles and the boundaries can produce shear 
bands 13, 12, 24\. Various theoretical models and simu- 
lations have pointed to the importance of the details of 
the competition between bubble-bubble dissipation and 
bubble-boundary dissipation [1, 0, [T3, 0] . A particularly 
interesting prediction exists for the case of the circular 
Couette geometry [J]. 

Circular Couette flow occurs for a material confined be- 
tween two concentric cylinders. One can drive flow from 
either the outer or inner boundary. An important char- 
acteristic of Couette flow (in the absence of any external 
dissipation) is the fact that the stress decays as 1/r 2 as 
one moves from the inner to the outer cylinder, indepen- 
dent of the driving. This establishes a natural situation 
for "classic" shear banding of a yield stress fluid, with 
flow occurring near the inner cylinder where the stress is 
the largest. An interesting feature of studies with bubble 
rafts is that this classic case is generally not observed. 
Instead, the shear bands are either discontinuous in the 
rate of strain or examples of discrete flow 0, 03 ■ I n this 
paper, we consider the case where external dissipation 
does exist, which modifies the spatial dependence of the 
internal stress, and explore the qualitative features of the 
shear banding. 

Recently, theoretical studies of the impact of dissipa- 
tion from an external boundary have predicted a range 
of striking behavior in Couette flow driven by the outer 
cylinder [4] . For sufficiently high external dissipation, one 
should observe a series of transitions from a single shear 
band at the inner cylinder to shear bands at the inner 
and outer cylinder, and ultimately for high enough ve- 
locities of the outer cylinder, a shear band exclusively at 
the outer cylinder. This effect is possible because of the 
velocity dependent dissipation of the confining bound- 
aries. In this paper, we report on an experimental con- 
firmation of the transition from single shear bands to two 
shear bands in the flow of a bubble raft. 

Our experimental setup utilizes bubble rafts in a Cou- 
ette geometry. We are able to run both with and without 
a top boundary for the bubbles. For the bubble raft, we 



find that the flow is driven predominately from the ro- 
tating outer boundary, and the drag due to viscous dissi- 
pation between the bubbles and the aqueous subphase is 
relatively negligible. For our Couette system, there is a 
circular dish (160 mm diameter, 20 mm deep) mounted 
from below on a shaft. A corrugated disk, approximately 
45 mm in diameter was set concentrically within the cir- 
cular dish. The shaft was supported by a set of bearings 
and connected using gears and driving bands to a vari- 
able speed stepper motor. The circular dish was filled 
with a homogeneous solution of high purity water, glyc- 
erol and miracle bubble solution in a 80:15:5 volume ra- 
tio. A stream of nitrogen gas was bubbled through the 
solution using a needle to generate a bubble raft in the 
spacing between the walls of the circular dish and inner 
disk. The bubbles in this 2D foam are polydisperse, with 
a diameters ranging between 2 mm to 9 mm. The cir- 
cular dish was rotated at different angular velocities in 
order to generate differing amounts of shear within the 
2D foam. The bubble raft was imaged from above using 
a standard CCD camera and digital images were saved to 
a computer every 0.1 s for a total of 4000 images. A stan- 
dard PIV technique was used to track the bubbles and 
compute the azimuthal velocity profile as a function of 
radial position. The velocity profiles were found to have 
converged over the observation period for all the angular 
velocities studied. An image of the apparatus is given in 
Fig. HI 

The experiment was run in two configurations: a) with 
the bubble raft open to the atmosphere and b) with a 
1/4 inch thick, transparent polycarbonate sheet in con- 
tact with the bubbles on top. With the addition of the 
top plate, a relatively uniform viscous drag is generated 
between the bubbles and the top plate. This drag force 
is velocity dependent. In both cases, the bubbles at the 
inner barrier are held stationary at the corrugated bar- 
rier. In case (a), slip at the outer barrier is eliminated by 
the insertion of 2 mm thick protrusions along the outer 
wall. In case (b) , the viscous drag between the top plate 
and the foam was stronger than that with the outer wall, 
and we were unable to eliminate slip between the outer 
wall and the last layer of foam. In this case, a corru- 
gated barrier was found to be insufficient to hold the 
bubbles in place, instead causing artifacts by frothing of 
the bubble raft even at low angular speeds. Therefore, 
in case (b) , the outer wall was moved at a higher angular 
velocity than the bubbles in contact with it. The vis- 
cous drag between the bubbles and the liquid subphase 
is not a significant contribution to the dynamics of the 
flow in comparison to that of the top plate and the outer 
wall. This is verified by checking that the bubbles do 
not undergo any rotational motion when only in contact 
with the inner barrier and the rotating subphase (without 
touching the outer wall, or the upper plate). A limitation 
in our experiment is that at high angular velocities with 
a top plate, a frothy solution was formed at the outer 
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FIG. 1: The top two panels are images of the bubble raft with 
(top panel) and without (middle panel) a top. The distance 
from the center of the trough to the outer wall is 16 cm. The 
lower panel is a photograph of the apparatus illustrating the 
experimental geometry. 



barrier that changed the dispersity of the foam from its 
initial condition. In the data presented in this paper, we 
limit our angular velocities to a regime where the froth 
formation does not occur significantly. The qualitative 
difference in the flow profile was not seen to be strongly 
influenced by the polydispersity of the bubble rafts used. 

The average velocity profile of the bubbles in our ex- 
periments have a circular symmetry reflecting the geom- 
etry of the shear. The average radial velocity is zero dur- 
ing steady shear as expected from the continuity equation 
in this radially confined system. We therefore focus on 
the azimuthal velocity (vg(r)) as a function of radial po- 
sition at different driving rates imposed. When rotating 
the outer boundary, a region of "no-flow" corresponds to 



rigid body rotation, vg(r) — Qr, where f2 is the angular 
rotation rate of the outer boundary. Therefore, to sim- 
plify the interpretation of the results, we present reduced 
velocities given by v(r) — v$(r)/ftr, and v(r) = constant 
is a region of "no flow". In Fig. [21 we plot v(r) for the 
case without a top. It should be noted that instead of 
using Q for the outer cylinder, we use f2 of the layer of 
bubbles adjacent to the outer cylinder. In the case of no 
top, these two numbers are equivalent because there is 
no slip. However, the presence of a top introduces slip, 
and for consistency we found it best to normalize by the 
first layer of bubbles. 

The measurements of v(r) in the absence of a confining 
top plate provide a comparison with previous measure- 
ments using a different Couette apparatus. The main 
feature of these measurements is the existence of a shear 
band at the inner cylinder for all rotation rates. However, 
it is worth pointing out that in this regime, the width of 
the the shear bands is independent of the rotation rate 
of the outer cylinder. 

If one defines critical radius for the shear band by the 
location at which the flow is indistinguishable from rigid 
body rotation, one observes a critical radius on the order 
of 10 of the smallest bubble diameter. Alternatively, one 
can fit the velocity profile in the shear band region to 
an exponential, and define a decay length as in Ref. jij. 
The results from this fit are show in Fig. 2J As with the 
critical radius, the decay lengths are independent of the 
angular velocity of the outer layer of bubbles. Previous 
studies of bubble rafts without a top have identified a 
flow regime for which the critical radius is independent 
of the rotation rate of the outer barrier This occurs 
for rates of strain below a critical value on the order of 
0.07 s _1 , and it is referred to as the discrete regime. Simi- 
lar behavior has been observed in three-dimensional foam 
systems [Uj . As the angular velocity of the outer cylin- 



der is on the order of the rate of strain, the observation 
of a rate of strain independent decay length is consistent 
with previous studies. 

With the addition of the top plate, the profiles at low 
rates of strain are similar to those without a top in that 
the flow is localized at the inner cylinder. However, if 
one considers the decay length for the shear band, one 
observes a dependence on the outer cylinder velocity. (As 
with the case without a top, an exponential fit was used 
to compute the decay length.) As the velocity increases, 
the decay length appears to saturate. This suggests that 
the dissipation from the top qualitatively changes the 
behavior of the system and plays an important role. For 
example, this behavior is completely different from what 
is observed in the discrete regime, and suggests a contin- 
uum model, such as the one reported on in Ref. Q , may 
be valid. Quantitative comparisons to the predictions re- 
ported in Ref. Q are severely limited, in part due to the 
fact that the magnitude of the dissipation between the 
plate and the bubbles is unknown, and that the model 
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FIG. 2: (color online)Plot of v(r) for bubbles subject to Cou- 
ette shear by driving the outer boundary at a fixed rotation 
rate. In this case, the bubbles are open to the atmosphere. 
The legend indicates the angular velocity (s _1 ) of the outer- 
most layer of bubbles during the applied flow. 



1 

§ o 

CO 

§ ° 

"8 o 

N 

ro 
E 

o 




4 5 
Radius (cm 



-■-0.0105 

-•-0.0188 
0.0281 

-T- 0.0377 
-J— —"-0.0836 
6 -► 0.0882 



FIG. 3: (color online) Plot of v(r) for bubbles subject to Cou- 
ette shear by driving the outer boundary at a fixed rotation 
rate. In this case, the bubbles are confined by a top plate. The 
legend indicates the angular velocity (s _1 ) of the outermost 
layer of bubbles during the applied flow. 



is relatively simplified. However, it is worth pointing out 
that the decay length for the case of rotating the inner 
barrier exhibits similar behavior in the model as observed 
in our experiments (see Fig. 8 of Ref. Q). However, this 
comparison is only indirect, as we rotate the outer bar- 
rier. 

As the velocity of the outer barrier is increased, we 
observe the behavior predicted in Ref. In addition 
to flow at the inner barrier, we observe flow at the outer 
wall. The region in between these two zones moves as 
a rotating rigid body. This results in a shear band at 
both boundaries. We were unable to sample higher shear 
rates due to limitations of our apparatus, and have not 
been able to observe the prediction of a single shear band 
at the outer wall. As with the shear bands at the inner 
cylinder, we fit the shear bands at the outer wall to an 
exponential profile and plot the decay lengths in Fig.Q] It 
is interesting to note that the trend for the outer barrier 
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FIG. 4: (color online) Plot of the decay length for the shear 
band as a function of the angular velocity (s~ ) of the out- 
ermost layer of bubbles during the applied flow. The results 
for the inner shear band are plotted versus the left-hand axis, 
with the circles for the system without a top and the squares 
for the system with a top. The results for the outer shear 
band with a top are plotted versus the right hand axis as 
triangles. 



decay lengths is opposite that of the one at the inner 
barrier. Though not reported on in Ref. [4], this is a 
useful test of the model that can be checked in detail in 
the future, and one can make a qualitative argument in 
the context of the model proposed in Ref. [4j . 

As suggested by Clancy et al. |4], the viscous damp- 
ing between the top plate and the bubbles is one source 
of the formation of shear bands. In the regime where 
this dominates, because the drag increases with increas- 
ing velocity, there is a tendency for the shear bands to 
move radially outwards with increasing shear rates. This 
observation is seen to hold at both the inner as well as 
the outer shear bands and is quantified by the increase 
in decay length for the inner shear band and decreasing 
decay band for the outer shear band with an increase in 
the rotation rate. It is also useful to note that while the 
viscous drag increases with radial position, the stress in- 
creases with decreasing radial position. The plateau be- 
tween the shear bands corresponds to a region where the 
viscous drag is balanced by the internal stresses within 
the bubble raft, resulting in a constant rotational veloc- 
ity. 

In summary, the observation of the dual shear bands 
represents agreement with the predictions of the theo- 
retical model described in Ref. [4j. For our system, the 
transition from a single shear band region at the inner 
boundary, to a dual shear band region at the inner and 
outer boundary occurs at a critical scaled azimuthal ve- 
locity for the outer row of bubbles of 0.03 s _1 . In the 
context of Ref. Q , this provide insight into the dissipa- 
tion due to the top plate. Our previous work has focused 
on two cases: bubble rafts with and without a top. For 
Couctte flow without a top, we have observed two regimes 
Q: discontinuous transitions and a discrete regime. As 
discussed, the results for the flow without a top reported 
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here are consistent with the discrete regime. As such, one 
does not expect the model discussed in Ref. 0] to be ap- 
plicable. Surprisingly, with a top, our current results are 
consistent with experiments from a range of groups that 
have confined bubbles either between two glass plates 
or between a fluid and a glass plate. In all of these 
cases, exponential profiles for the velocity are observed 
13 , 53 . 2a |. The differences between these studies is the 
behavior of the width of the shear band as a function of 
the applied rotation rate. However, these differences in 
rate dependence can be explained in the context of the 
continuum model (sec for example. Fig. 8 of Ref. 4] and 
Ref. [13] which extends the model of 0]). 

Finally, the direct observation of a transition from a 
single shear band to a double shear band as a function 
of increasing velocity of bubbles for a bubble raft in a 
Coucttc geometry confirms the prediction of the model 
described in Ref. [J] that such behavior occurs in the 
presence of sufficient external drag on a foam system. Ex- 
perimentally, the fact that this shear band occurs when a 
top plate is added confirms that confining boundaries in 
quasi-two dimensional foam systems provide substantial 
drag that must be understood and accounted for in the- 
oretical models. Though it remains clear from other ex- 
periments that this mechanism will not explain all cases 
of shear-banding, these experiments provide strong ev- 
idence that it is an important mechanism in situations 
involving external drag. Some interesting future direc- 
tions will be to look at the impact of polydispersity (as 
in Ref. [2l|]) an d to explore in more detail the impact 
of velocity dependent drag forces in the regime in which 
discontinuous shear localization exists. 
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